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The block operator matrix theory is used to investigate the problem of a single qubit. We will establish a 
connection between the Riccati operator equation and the possibility of obtaining an exact reduced dynamics 
for the qubit in question. The model of the half spin particle in the rotating magnetic field coupling with the 
external environment is discussed. We show that the model defined in such a way can be reduced to a time 
independent problem. 
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I. INTRODUCTION 

Exactly solvable models for decoherence play an im- 
portant role both in the theory of the open quantum 
system and quantum information theory^'-'. Unfortu- 
nately, most of the models describing the process of de- 
coherence can not be solved exactly. However, there is a 
wide class of models for which exact reduced dynamics ' 
is known. Those models deal with the case where the 
energy transfer between the system and the environment 
is not present. This phenomenon is known as pure deco- 
herence or dephasing'. It has been found that general- 
ization of the dephasing models to the case where energy 
is exchanged between the system and the environment is 
straightforward, but for most of this generalization ana- 
lytical solutions were not obtained. It is only natural to 
wonder why the dephasing models can be solved easily 
whereas even the most basic generalizations pose such a 
difficult task for the scientists. 

In this manuscript we show that obtaining the exact 
reduced dynamics of the model of one qubit interacting 
with the environment is at least as difficult as solving the 
Riccati operator equation associated with the Hamilto- 
nian defining the model. First we will discuss the pro- 
cedure allowing one to obtain the density matrix for the 
system using the block operator matrix perspective. 

The general form of the Hamiltonian describing the 
qubit Q coupling with the external environment (heat 
bath) E can be written as follows ' 



Hqe = HQ(g)lE + Iq'E> He + H,, 



(1) 



where Hq, He are the Hamiltionian of the qubit and 
the environment respectively and Hint represents the in- 
teraction between the systems. Hamiltonian Hqe acts 
on T-Lq ^He space, where 'Hq and He are the Hilbert 
spaces for the system and the environment respectively. 
For most models it is assumed that the initial state of 
the Q + E system has the following form pqe = pq® pe, 
which means that there is no correlation between Q and 



E initially (see ' and ref. therein). Our analysis is free of 
this assumption. The state of the Q system at any given 
time t takes the form: 



PQit)=TTEiUtPQ®PEU}), 



(2) 



where Ut is the evolution operator of the Q + E system 
and by Tr^(-) we denote the partial trace. The PQ{t) is 
called reduced dynamics (with respect to the degree of 
freedom of the environment). From now on the quantity 
PQ{t) will be called the solution of the model. In the 
case of Hq = and He = H, where H is the arbitrar- 
ily separable Hilbert space (in general dimH = oo) the 
following isomorphism holds C"^ ^ H = H (B H. There- 
fore, any given operator A acting on the ®H space 
can be thought of as the 2x2 block operator matrix 
(BOM) [Aij], where A^^, {i,j = 1,2) act on H. In this 
notation the procedure of calculating partial trace Tr e is 
very intuitive, namely 



TrEiA) 



TrAn Tr^is 
TrAai TrA22 



(3) 



^'Electronic mail: bartck.gardas@gmail.com 



where Tr(-) is a trace on H. One can easily see that 
obtaining reduced dynamics pQ(t) is very simple. How- 
ever, the equation (2) is far less useful than its theoretical 
simplicity might indicate. The reason is that one can not 
determine the exact block operator 2x2 matrix form 
of the evolution operator Ut of the system Q + E. The 
task becomes even more difficult when the Hamiltonian 
is time dependent. 

There have been few theories resolving the problem of 
finding reduced dynamics both for time dependent and 
time independent Hamiltonians''*'. Furthermore a ma- 
jority of the scientists focus their effort on a numerical 
methods and on perfecting the approximation methods ' . 
As a consequence most of the research on the quantum 
information theory is based on a numerical rather then 
an analytical approach. As a result during past several 
years no progress has been made in solving the known 
models. 

The main purpose of this manuscript is to present an 
analytical approach. We consider one of the most estab- 
lished and useful models, namely the spin 1/2 (qubit) in 
the rotating magnetic field. In the case where no coupling 
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with the external environment is present, an analytical 
solution can be found in an elegant and simple manner"^. 
If the mentioned coupling (modelled by quantum system 
of infinite number of degree of freedom) is present; how- 
ever, the exact solution has not been found yet. We will 
not address this in the current manuscript; however, we 
will show that this model can be effectively reduced to 
the time independent problem (Section II). Moreover, 
we will show that the solution of any given model with a 
time independent Hamiltonian requires solving the Ric- 
cati operator equation associated with the Hamiltonian 
H defining the problem (Section III). In other words, 
we will establish the connection between the problem of 
decoherence in physics and the mathematical problem of 
resolving the Riccati operator equation. Furthermore, 
using the results of Section III we will discuss the pos- 
sibility of obtaining an exact solution to the analyzed 
problem from the set of differential equations on "H "H 
(Section IV). Finally, in Section V we consider a spin- 
bozon model as an example. Section VI is a summary of 
the paper. 

II. SPIN HALF IN A ROTATING MAGNETIC FIELD 
AND IN CONTACT WITH ENVIRONMENT 

Let us consider a single qubit in rotating magnetic field 
interacting with its environment. The qubit-environment 
time-dependent Hamiltonian reads 

H{t,(3)=HQ[t,p)®lE + h®HE + H,^u (4) 

where Hglt, (3) and He are Hamiltonians of qubit Q and 
the environment respectively and Hint represents the in- 
teraction between Q and the environment. It is assumed 
that Hint takes the form /((T3) (E) V, where y is a Hcr- 
mitian operator acting on He and /(era) is an analytic 
fmiction of (T3. Hamiltonian HQ{t,j3) is given by 

Hgit, (3) = /3(T3 + a ((Ti cos [uot) + (T2 sin {ujt)) , (5) 

and it represents a spin system in rotating magnetic field 
B{t), where 

B{t) = [Bi cos (wO , Bi sin (wi) , Bq] . (6) 

Here, a ~ ^ Bi and [3 ~ ^cjo ^ Bq, where Bq, Bi 
are amplitudes of the magnetic field* ' . 

The model described by the Hamiltonian (4) cannot be 
solved exactly in this general case. By this we mean that 
the exact reduced dynamics pg (t) for that model are not 
known. Let us now focus on another model defined by 
the Hamiltonian H{I3) = H{0,13), where H{t,(3) is given 

by (4). 

We will show that if rjt is a solution to the model de- 
scribed by the Hamiltonian (4), and pt(/3) represents a 
solution of the model with Hamiltonian H{f3), then the 
following equation holds 

Vt = Vtpt{l3-'^)v,\ (7) 



where 

Ft = diag (e-*"*/2^ e''"*/2^ . (8) 

From equation (7) and (8) we see that if reduced dy- 
namics /9t(/3) is known then all one needs to do to obtain 
the solution to the model of the H{t, (3) Hamiltonian is to 
introduce an effective parameter /3e// := — ^, replace (3 
by /3e//j and perform a unitary transformation (8). Since 
the procedure explained above is very simple we can ef- 
fectively reduce the problem of solving model (4) to one 
of solving the model H{/3). 

In order to prove the equation (7) let us note that the 
Hamiltonian (4) satisfies the following condition (h ~ 1) 

H{t,(3)^e'^'H{P)e-'^\ (9) 

where K = —^a^ ® 1e- This can be easily proven using 
the Baker-Campbell-Hausdorff formula "'. As was shown, 
in^ every quantum system with Hamiltonian H{t, (3) sat- 
isfying (9) for some Hermitian operator K there evolves 

UtiP) = eK^e-^^^sm\ H,ff{l3) := H{P) + K. (10) 

Note that in general [H{f3),K] 7^ and therefore 
[Heff{l3),K] 7^ 0. In our case, from equation (4) we 
learn that H{(3) = {Pa^ + aai) ^ 1e, thus 

Heffif3) = iPa3 + aai)®lE- '^(73 ®1e (H) 
= ((/?- I) era + 



From equations (10) and (12) we have 

[/,(/3)=e*^*[/t(/3-f), (12) 

where Ut{(3) is the evolution operator generated by 
H{t,(3). Let pt{l3) and fit be a density operator for the 
closed system Q + E associated with Hamiltonian i?(/3) 
and H(t, (3) respectively in arbitrary time t. Let us also 
assume that po{l3) ~ fio = p. Using equation (12) one 
can easily see that 

f,t = Ut{f3)pU}{l3) (13) 
= e*^*f/t(/3 - f )pf//(/3 - f )e-*^* 
= VtPt{l3-^)V,\ 

where we introduced Vt = e'-^*. To end the proof we 
wiU show that if ^1, A2 G B{'H®'H) are a 2 x 2 block 
operator matrix of the form Ai = (g) 1^, (i = 1, 2) and 
B = [By ] G Bin © n) then 

TrE(iiBi2) = AiTve{B)A2. (14) 

Equation (14) follows from the linearity of trace Tr op- 
eration and definition (3) of partial trace. Note that 
Vt = Vt®lE, where Vt is given by equation (8), thus 
taking partial trace of equation (13) and using (14) we 
obtain (7) with Vt given by (8). 
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III. OPERATOR RICCATI EQUATION 

So far we have shown that the solution rjt can be easily 
constructed from ptiP). Now, we will pay attention to 
the possibility of obtaining an exact solution pt{P)- Let 
us now rewrite Hamiltonian H{j3) as a block operator 
matrix^ ^ 



H 



+ , /3 a 



(15) 



where we introduced H± = He ± V. Since Hamilto- 
nian (15) is time-independent, we can write the evolution 
operator as Ut = exp {—iH{(3)t). We see that the main 
problem here is how to write down [/j as 2 x 2 BOM. 

If a = this problem is trivial. On the other hand for 
a the diagonalization of 2 x 2 BOM is required which 
is not a trivial problem^-. With every Hermitian 2x2 
BOM of the form: 



R 



A B 
St C 



A,B,C<En 



(16) 



we can associate the operator Riccati equation^ '' 

XBX + XA- CX - =0, (17) 

where X E H. Solution X of the equation (17), if it 
exists can be used to construct 2x2 BOM: 



in such a way 

U^'RUx 



Ie -X^ 
X Ie 



A + BX 
C-Bt^t 



(18) 



(19) 



From the above consideration we see that to diagonalize 
Hamiltonian (15) we have to solve the following Riccati 
equation: 



aX^ + X{H+ +13)- {H_ - f3)X 



0. 



(20) 



Unfortunately, we do not know how to do that. Note 
that if a = then X = is a solution. This is obvious 
since in that case H{/3) is already in the diagonal form. 
Note also that even if /3 = this problem is still very 
complicated. 



that we circled back to the point where we started since 
in writing the state as a column vector we need to 
diagonalize the matrix H{/3). Nothing could be further 
from the truth. To see this let us introduce operators U 
and Jt in the following way: 



U- 



1 

71 



1 1 

1 -1 



Jt = cxp^iaa^t) . 



(21) 



Let us also define |4't) = JtU\'^t), and we can easily see 

that i\'^t) = Ht\^t), where the periodic Hamiltonian is 
given by 



He z;{V + I3) 
zt{V + l3) He 



Zt 



(22) 



The Riccati equation associated with Ht reads 

Xiz;Vp)X + XHe - HeX - ztVp = 0, (23) 

where Vp = V + /3. Straightforward calculations show 
that Xt = Zt is a solution of the Riccati equation (23). 
According to (19) we have 



SjHtSt = 



H++(3 
-/? 



(24) 



where St = '^Uzt and Uz^ is given by (18), namely: 



1 

Zt 



1 



(25) 



Note that St is a unitary 2x2 block operator matrix. 
We see that if one could solve the Schrodinger equation 
for I'J't) then our problem would be solved. Formally, we 
can always do that using chronological operator T, the 
solution is given by^^ 



l^-t) =Tcxp(-Z / HrdT)\^'o) 




(26) 



Sadly, the form (26) of the solution has little use due to 
the presence of the chronological operator T. Neverthe- 
less, it is interesting to note that the connection between 
the models (15) and (24) is well define and the solution 
to the equation (23) can be easily found, yet finding the 
solution to the equation (20) poses a big problem. 



IV. DIFFERENTIAL EQUATION APPROACH 

Let us now transform the problem of solving a Riccati 
equation (20) to the problem of solving a Schrodinger 
equation on 'H © 7^, with the Hamiltonian given by (15). 
Let l^-t) = [|V't),|0t)]*, then l^-t) satisfy zj^-t) = 
H{j3)\'it)- Of course, we can always write = 
exp(— ii?(/3)t)|\E'o), but this form of the solution is useless 
since Ut does not have a 2 x 2 BOM form. It may seem 



V. EXAMPLES 

Up until now we did not choose the specific form of the 
operators He and V, which means that the analysis we 
presented was very general. That fact implies an impor- 
tant concept, namely that the analysis that we carried 
out does not depend on the particular choice of a heat 
bath. It is crucial, however that the coupling of the qubit 
with the environment is given by the matrix /((Ts), where 
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/ is an arbitrary analytical function. It is interesting to 
consider the model where operators He and V = V{g) 
are defined as follows: 



He = / du) (jj a' {uj)a{u}), 



(27) 



where 0^(0;) and a{ui) are bozonic annihilation and cre- 
ation operators respectively and they satisfy the commu- 
tation relations: [a(w), a''^(a;')] = S{uj — oj'), uj, uj' > 0. 
V{g) is given by 



If X = xIe, where x G C, then one can write the above 
equation as 

(mi2X^ + (mil - 77122)2; - 771^2) V = 0, (35) 
or equivalently as 



mi2x'^ + (77711 — m22)x — ~ 0. 



(36) 



As a result, we see that for the dephasing case the Riccati 
equation simplifies to the quadratic equation and there- 
fore solution X can be easily found. 



V{g) ^ dw {g*{uj)a{uj) + g{uj)a^uj)) , 



(28) 



where g € L^[0,oo]. Operators He and V{g) given 
by (27) and (28) define the bozonic heat bath ' of the 
qubit. One can find that 

H^ = W{g)HEW{gy +Cig), (29) 
H_^Wig)^HEWig)+Cig), (30) 

where C{g) is a certain constant. One can always rescale 
the Hamiltonian so that C{g) — 0, thus we will omit 
the constant. The unitary Weyl's operator has the form 
W{g) = cxp{A{g)), where 



A{g) = / {g*{u;)a{uj) - g{uj)a^uj)) 



(31) 



In the case of the a ~ model can be solved exactly^ . If 
a 7^ obtaining the exact reduced dynamics, according 
to (20) is at least as difficult as solving the equation (to 
simplify we put (3 = 0) 



aX^ + XiWHEWf) - {W^iHeW^X ■ 



0. (32) 



The solution of the equation (32) is yet to be discovered. 

As a second examples let us consider a pure decoher- 
ence case. In this situation [Hq (g) lE,Hint] = 0. Let 
Hint = M ^V, where M is a arbitrary Hermitian 2x2 
matrix. Since operators Hq ^ 1e and M (S)V commute, 
we need to diagonalize the following matrix to solve our 
problem: 



I2(E)He + M(E)V 



HE + rriiiV 77712V 

"^12^ He + 7)7.22 V" 



(33) 



The Riccati equation associated with BOM (33) takes 
the form 

mi2XVX + X{HE + mi2V) (34) 

- (i?B+ 777 2 2^)^- 777 * 2^ = 0. 



VI. SUMMARY 

111 this paper the problem of the exact solution of the 
decoherence model has been connected to the Riccati op- 
erator equation. It was shown that obtaining the exact 
reduced dynamics is as problematic as resolving the Ric- 
cati equation to say the least. Furthermore, we simplified 
a wide class of problems described by the time dependent 
Hamiltonian to the time independent problems. One can 
easily learn from this paper that solving the time depen- 
dent Riccati (23) equation is very simple. 

We strongly believe that solving the model we analyzed 
is crucial and that it can contribute to the progress and 
verification of the adiabatic theorem for open quantum 
systems ' in analogy to the contribution of the half spin 
particle model with Hamiltonian Hq (i, /?) to the progress 
of the adiabatic theorem for the quantum closed systems. 
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